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Using (223.7 + 1.4) x 10° J/w events accumulated with the BESIII detector, we study rj decays 
to ġġ and wọ final states. The branching fraction of re — ¢¢ is measured to be Br(ne — $$) = 
(2.5 + 0.3*03 + 0.6) x 107°, where the first uncertainty is statistical, the second is systematic, and 
the third is from the uncertainty of Br(J/w — yne). No significant signal for the double Okubo- 
Zweig-lizuka suppressed decay of 7. + w is observed, and the upper limit on the branching fraction 
is determined to be Br(n. > wọ) < 2.5 x 1074 at the 90% confidence level. 


PACS numbers: 13.25.Gv, 13.20.Gd 


I. INTRODUCTION 


Our knowledge of the 7. properties is still relatively 
poor, although it has been established for more than thir- 
ty years |l]. Until now, the exclusively measured decays 
only sum up to about 6396 of its total decay width [2]. 
The branching fraction of ne — $9 was measured for the 
first time by the MarkIII collaboration and improved 
measurements were performed at BESII [5 with a pre- 
cision of about 40%. The decay ne — wọ, which is a dou- 
bly Okubo-Zweig-lizuka (OZI) suppressed process, has 
not been observed yet. 

Decays of r) into vector meson pairs have stood as a be- 
wildering puzzle in charmonium physics for a long time. 
This kind of decay is highly suppressed at leading order 
in Quantum Chromodynamics (QCD), due to the helic- 
ity selection rule (HSR) ie. Under HSR, the branching 
fraction for ne — $ó was calculated to be ~ 2 x 1077 R. 
To avoid the manifestation of HSR in charmonium de- 
cays, a HSR evasion scenario was proposed [8]. Improved 
calculations with next-to-leading order I and relativis- 
tic corrections in QCD yield branching fractions varying 
from 107? to 1074 mi Some nonperturbative mech- 
anisms, such as the light quark mass corrections [12], the 
3 P) quark pair creation mechanism and long-distance 
intermediate meson loop effects [14], have also been phe- 
nomenologically investigated. 

However, the measured branching fraction, Br(r. > 
oo) = (1.76 + 0.20) x 107? [A [15], is much larger than 


those of theoretical predictions. To help understand 
the 7. decay mechanism, high precision measurements 
of the branching fraction are desirable. In this paper, 
we present an improved measurement of the branching 
fraction of ne — ġġ, and a search for the doubly OZI 
suppressed decay ». — wo. The analyses are performed 
based on (223.7 + 1.4) x 106 J/w events [16] collected 
with the BESIII detector. 


II. DETECTOR AND MONTE CARLO 
SIMULATION 


The BESIII experiment at BEPCII is an upgrade 
of BESII/BEPC . The detector is designed to study 
physics in the 7-charm energy region [19]. The cylindri- 
cal BESIII detector is composed of a helium gas-based 
main drift chamber (MDC), a time-of-flight (TOF) sys- 
tem, a CsI (TI) electromagnetic calorimeter (EMC) and a 
resistive-plate-chamber-based muon identifier with a su- 
perconducting magnet that provides a 1.0 T magnetic 
field. The nominal geometrical acceptance of the detec- 
tor is 93% of 4» solid angle. The MDC measures the 
momentum of charged particles with a resolution of 0.5% 
at 1 GeV/c, and provides energy loss (dE/dx) measure- 
ments with a resolution better than 696 for electrons from 
Bhabha scattering. The EMC detects photons with a res- 
olution of 2.5% (5%) at an energy of 1 GeV in the barrel 
(end cap) region. 
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To optimize event selection criteria and to understand 
backgrounds, a GEANT4-based Monte Carlo (MC) 
simulation package, BOOST, which includes the descrip- 
tion of the geometries and material as well as the BESIII 
detection components, is used to generate MC samples. 
An inclusive J/~-decay MC sample is generated to study 
the potential backgrounds. The production of the J/w 
resonance is simulated with the MC event generator 
KKMC , while J/w decays are simulated with BESEVT- 
GEN [22] for known decay modes by setting the branching 
fractions to the world average values [2], and with LUND- 
CHARM for the remaining unknown decays. The anal- 
ysis is performed in the framework of the BESIII offline 
software system [24]. which handles the detector calibra- 
tion, event reconstruction and data storage. 


III. EVENT SELECTION 


The 7. candidates studied in this analysis are produced 
by J/w radiative transitions. We search for ne > $$ and 
wọ from the decays J/ > yọġ and ywd, with final 
states of y2(K* K ) and 3yK * K m^. , respectively. 
The candidate events are required to have four charged 
tracks with a net charge of 0, and at least one or three 
photons, respectively. 

Charged tracks in the polar angle region | cos 0| < 0.93 
are reconstructed from the MDC hits. They must have 
the point of closest approach to the interaction point 
within +10 cm along the beam direction and 1 cm in the 
plane perpendicular to the beam direction. For the par- 
ticle identification (PID), the ionization energy deposited 
(dE/dx) in the MDC and the TOF information are com- 
bined to determine confidence levels (C.L.) for the pion 
and kaon hypotheses, and each track is assigned to the 
particle type with the highest PID C.L. For the decay 
J/i — qw > 9y K * K mm, two identified kaons are 
required within the momentum range of 0.3—0.9 GeV 
with an average efficiency of about 896. For the decay 
J/ > yp > 42(K TK), no PID is required. The in- 
termediate states, @ and w, are selected using invariant 
mass requirements. 

Photon candidates are reconstructed by clustering en- 
ergy deposits in the EMC crystals. The energy deposited 
in the nearby TOF counters is included to improve the 
photon reconstruction efficiency and energy resolution. 
'The photon candidates are required to be in the bar- 
rel region (|cos0| < 0.8) of the EMC with at least 25 
MeV total energy deposition, or in the end cap regions 
(0.86 « | cos0| « 0.92) with at least 50 MeV total energy 
deposition, where 0 is the polar angle of the photon. The 
photon candidates are, furthermore, required to be sepa- 
rated from all charged tracks by an angle larger than 10? 
to suppress photons radiated from charged particles. The 
photons in the regions between the barrel and end caps 
are poorly measured and, therefore, excluded. Timing 
information from the EMC is used to suppress electronic 
noise and showers that are unrelated to the event. 


Kinematic fits, constrained by the total e*e^ beam 
energy momentum, are performed under the J/ — 
y2(K*K-) and 3yK* K m^7- hypotheses. Fits are 
done with all photon combinations together with the four 
charged tracks. Only the combination with the smallest 
kinematic fit xç is retained for further analysis, and 
Xác < 100 (40) for J/j > 42(K* K ) (24K* K mm) 
is required. These requirements are determined from MC 
simulations by optimizing S/V S + B, where S and B are 
the numbers of signal and background events, respective- 
ly. 

Two ó candidates in the J/v — ydd decay are re- 
constructed from the selected 2(K* K^) tracks. Only 
the combination with a minimum of IME p- — Map? + 
IMO). — MsJ? is retained, where M, (i = 1,2) 
and M, denote the invariant mass of the K* K^ pair 
and the nominal mass of the ¢-meson, respectively. A 


scatter plot of M I g- versus M pel g- for the surviving 
events is shown in Fig. [] (a). There is a cluster of events 
in the $$ region indicated as a box in Fig. [I] (a) orig- 
inating from the decay J/W — y¢¢. Two ¢ candidates 
are selected by requiring |Mg+g- — M| < 0.02 GeV/c”, 
which is determined by optimizing S/VS + B, also. 

For the decay J/ — ywo — yK* K-n*n T7, the 
photon combination with mass closest to the 7° nominal 
mass is chosen, and |M; — M;o| < 0.02 GeV/c? is re- 
quired. A scatter plot of the Mg g- versus Mj, 4o 
for the surviving events is shown in Fig. [I] (b). Three 
vertical bands, as indicated in the plot, correspond to 
the m, w and @ decays into m+a~7°, and the horizon- 
tal band corresponds to the decay ¢ + K* K-. For the 
selection of J/4h — ywọ candidates, the ó and w require- 
ments are determined, by optimizing S/V S + B, to be 
|Mr+r-70 — Mwl < 0.03 GeV/c? and |Mi«y- — Mg| < 
0.008 GeV/c?. 


IV. DATA ANALYSIS 
A. Observation of n. > %% 


Figure 2] shows the invariant mass distribution of the 
¢¢-system within the range from 2.7 to 3.1 GeV/c?. The 
Te Signal is clearly observed. Background events from 
J/w decays are studied using the inclusive MC sample. 
The dominant backgrounds are from the decays J/w > 
yo6K* K- and J/y > 4K* K K*K- with or without 
an rj, intermediate state, which have exactly the same fi- 
nal state as the signal, and are the peaking and non peak- 
ing backgrounds in the 2(K* K ^) invariant mass distri- 
bution. In addition, there are 43 background events from 
the decays J/w > #f1(1420)/f1(1285) with fı decay to 
KEK’ and J/v — oK*(892)* K* with K*(892)* de- 
cay to Kz) which have a final state of 3?2(K * K) 
similar to that of the signal. These background decay 
channels have low detection efficiency (« 0.196), and do 
not produce a peak in the T< signal range. The expected 
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Fig. 2: Projection of fit results onto the M; spectrum. The 
dots with error bars denote the data, the solid line histogram 
is the overall result, the dot-dashed histogram is the m. sig- 
nal, the filled red histogram is the combined backgrounds es- 
timated with exclusive MC simulations, the dotted histogram 
denotes non 75. decays, and the long-dash histogram is the 
interference between the ne and non ne decays. 


yields of background events are 26 and 75 for the peaking 
and non peaking backgrounds, respectively, determined 
with MC simulation. As a cross-check, the backgrounds 
are also estimated with the events in the ¢ sidebands re- 
gion in data, and then using the MC information of the 
Ne > K+K- and 2(K* K ) to scale the rj. events in 
boxes B, C and D to the signal region A, and total 104 
events are obtained. 


To determine the 7. — ¢@ yield, an amplitude analy- 
sis is performed on the selected candidate 1,276 events. 
We assume the observed candidates are from the process 
J/ — yo with or without the ry. intermediate state in 
the ġġ system. The amplitude formulas are construct- 
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for the decay J/y — y2(Kt K7), and (b) Mg+g- versus M, o for the 


ed with the helicity-covariant method [23], and shown in 
the appendix. The 7. resonance is parametrized with the 
Breit- Wigner function multiplied by a damping factor 


1 Z(E) 


Ms) = m. iMT FBS)’ 


where s is the square of ġġ invariant mass, and M and 
I are the 7. mass and width, respectively. The damping 
2 


factor is taken as F(E,) = exp(- 322) with 8 = 0.065 
GeV lea, and the photon energy EP corresponds to the 
vs — M. 

In the analysis, the decay J/w — yne > yoo and the 
nonresonant decays J/w — yd with different quantum 
numbers J^ (spin parity) in the $$ system are taken 


into consideration. The differential cross section do /dQ 
is calculated with 


da 
ADU > [An (Ao, Ay, Ài, A2) 


helicities (2) 
+ Y Ago Ay At, AP, 
JP 


where A,,, is the amplitude for the J/(Ao) > (Ay) Ne > 
7(A1)¢(A2), with the joint helicity angle Q, and AT, is 
the amplitude for the nonresonant decay J/w —> yoo 
with JP for the ġġ system. To simplify the fit, only 
the nonresonant components with J^ = 0+, 07 and 2* 
are included, and the components with higher spin are 
ignored. The symmetry of the identical particles for the 
ó$-meson pair is implemented in the amplitude. 

The magnitudes and phases of the coupling constants 
are determined with an unbinned maximum likelihood 
fit to the selected candidates. The likelihood function 
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for observing the N events in the data sample is 


N 


£=]] Pes, (3) 


i=l 


where P(z;) is the probability to observe event ¿ with 
four momenta x; = (py,p4,P¢)i, Which is the normalized 
differential cross section taking into account the detection 
efficiency (€i), and calculated by 


(do /dQ);€; 


Finge OMC 


(4) 
where the normalization factor omc can be calculated 
by a signal MC sample J/u — yé¢ with Nuc accepted 
events. These events are generated with a phase space 
model and then subjected to the detector simulation, and 
passed through the same events selection criteria as ap- 
plied to the data. With a MC sample which is sufficiently 
large, omc is evaluated with 


1 a (5) 
Nuc i=l dQ "i 


For a given N events data sample, the product of e; 
in Eq.(3) is constant, and can be neglected in the fit. 
Rather than maximizing £, 7 = — ln £ is minimized us- 
ing MINUIT [27]. 

In the analysis, the background contribution to the 
log-likelihood value (In Lbkg) is subtracted from the log- 
likelihood value of data (In Laata), ie. In £ = InLaata — 
In Lbkg, where In Lbkg is estimated with the MC simulat- 
ed background events, normalized to 101 events including 
peaking and nonpeaking 7. background. 

In the fit, the mass and width of rj. are fixed to the pre- 
vious BESIII measurements Ps, Le. M = 2.984 GeV/c? 
and T = 0.032 GeV. The mass resolution of the ry. is not 
considered in the nominal fit, and its effect is considered 
as a systematic uncertainty. The fit results are shown in 
Fig. | where the rightmost peak is due to backgrounds 
from J/w > K+K- decay. The m yield from the fit is 
Nn. = 549 + 65, which is derived from numerical integra- 
tion of the resultant amplitudes, and the statistical error 
is derived from the covariance matrix obtained from the 
fit. 

To determine the goodness of fit, a global x: is calcu- 
lated by comparing data and fit projection histograms, 
defined as 


TMC = 


N (NPE _ NES 


5 
Xe = DIG, with xj = > — 
j=1 


3 6 

i=1 N. A. s 
where NRT and Nj;** are the numbers of events in the ith 
bin of the ith kinematic variable distribution. If N a is 
sufficiently large, the X2 is expected to statistically follow 
the x? distribution function with the number of degrees 
of freedom (ndf), which is the total number of bins in 


histograms minus the number of free parameters in the 
fit. In a histogram, bins with less than ten events are 
merged with the nearby bins. The individual x7 give a 
qualitative evaluation of the fit quality for each kinematic 
variable, as described in the following. 

Five independent variables are necessary to describe 
the three-body decay J/« — yp. These are chosen to 
be the mass of the ¢¢-system (Mo), the mass of the y¢- 
system (Mo), the polar angle of the y (04), the polar an- 
gle (05) and azimuthal angle (yy) of the ¢-meson, where 
the angles are defined in the J/w rest frame. Figure 
shows the comparison of the distributions of My and 
angles between the global fit and the data. A sum of 
all of the x7 values gives x2 = 215 with ndf=191. The 
quality of the global fit (x2/ndf) is 1.1, which indicates 
good agreement between data and the fit results. 

'To validate the robustness of the fit procedure, a pseu- 
dodata sample is generated with the amplitude model 
with all parameters fixed to the fit results. A total of 
2936 events are selected with the same selection criteria 
as applied to the data. An identical fit process is car- 
ried out, and the ratio of output ne signal yield to input 
number of events is 1.03 + 0.03. 


B. Search for n. > wọ 


Figure [4] shows the w@ invariant mass distribution in 
the range from 2.70 to 3.05 GeV/c? for the selected can- 
didate events of J/4» — ywd, and no significant ne signal 
is observed. The background events from J/w decays are 
dominated by J/w — nó with n” > yw. A small amount 
of background is from the decays J/v — fo(980)w > 
K*K-w and J/ — fxw — n? K*K-w, where fx 
stands for the f1(1285) and /1(1420) resonances. The 
sum of all above backgrounds estimated from inclusive 
MC samples is small compared to the total number of 
selected candidates and appears as a flat Mwg distribu- 
tion, as shown in Fig. [4] 

To set an upper limit for the branching fraction 
Br(ne — wọ), the signal yield is calculated at the 90% 
C.L. by a Bayesian method |2], according to the distribu- 
tion of normalized likelihood values versus signal yield, 
which is obtained from the fits by fixing the 7, signal 
yield at different values. 

In the fit, the shape for the rj. signal is described by the 
MC simulated line shape by setting the mass and width 
of n. to the BESIII measurement ; the known back- 
ground estimated with MC simulation is fixed in shape 
and magnitude in the fit; and the others are described by 
a second-order Chebychev function with floating param- 
eters. The distribution of normalized likelihood values is 
shown in Fig. B] and the upper limit of signal yield at 
the 9096 C.L. is calculated to be 18. 

To check the robustness of the event selection crite- 
ria, especially the dependence on Br(n. — wọ), the re- 
quirements of kinematic fit Y? and ¢/w mass windows 
are reoptimized with the measured upper limit. The ne 


EVENTS / 10 MeV/c2 


EVENTS / 0.1 


EVENTS / 0.1 


EVENTS / 0.1 


120 


100 


c 


P 


Fig. 3: Distributions of (a) the yọ invariant mass M.,5; (b) the polar angular of the photon cos 84; (c) the polar angular of $ 
mesons cos 05; (d)and the azimuthal angular of ó mesons yy. The dots with error bar are the data, the solid line histograms 
represent the total fit results, and the filled histograms are the non-J/1» — y¢¢ backgrounds estimated with the exclusive MC 
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Fig. 4: Results of the best fit to the Mwg distribution. Dots 
with error bars are data, the solid curve is the best fit result, 


corresponding to a 7c signal yield of 104 


- 6 events, the shaded 


histogram is the background estimated from exclusive MC 
samples, the dashed curve indicates the 7. signal, and the 


dotted curve is the fitted background. 
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Fig. 5: Normalized likelihood distribution versus the 7. yield 
for nc > wọ. 


V. 


SYSTEMATIC UNCERTAINTIES 


The following sources of systematic uncertainties are 


considered in the measurements of branching fractions. 


1. Number of J/w events 
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The number of J/z events is determined using its 
hadronic decays. The uncertainty is 0.6% [16]. 


. Photon detection efficiency 


The soft and hard photon detection efficiencies are 
studied using the control samples Y’ — z9z9J/40, 
with J/4» decay e*e- or u* i^. and J/ — pr > 
v^ 7, respectively. The difference in the photon 
detection efficiency between the MC simulation and 
data is 196, which is taken as the systematic uncer- 
tainty. 


. Kaon/pion tracking and PID efficiency 


The uncertainties of kaon/pion tracking and PID 
efficiency are studied using the control samples 
J/i > m! pp and J/y > KZ Kr, with the decay 
K9 — m m+ The uncertainties for tracking 
and PID efficiencies are both determined to be 196 
per track. 


. Branching fractions 


The uncertainties of branching fractions for J/w — 
Me, $ — K* K- , and wu > ntron’ are taken from 
the PDG [2]. 


. Kinematic fit 


To estimate the uncertainty associated with the 
X? requirement of the kinematic fit for the final 
state y2(K*K—), we select the candidate events 
of J/i — yọ¢ by requiring x? < 20, 60 or 150, 
and the T signal yields are reestimated with ampli- 
tude analysis. The largest deviation to the nominal 
branching fraction, 6.796, is taken as the systematic 
uncertainty. 


For the final states y K * K ntr 7°, we redeter- 
mine the upper limit on the branching fraction 
with the alternative requirement of the kinematic 
fit x? < 20,30,50 or 60, and the largest deviation 
to the nominal value, 2.496 at x? « 30, is taken as 
the systematic uncertainty. 


. Mass window 


The uncertainties associated with the ¢/w mass- 
window requirement arise if the mass resolution is 
not consistent between the data and MC simula- 
tion. The uncertainty related to the $-mass window 
requirement is determined with the control sample 
y > "yXeJ: XeJ > oo, and $ — K*K-. The 
difference in $-selection efficiency is estimated to 
be 0.7% and 1.1% for the ne — $9 and ne > wọ 
modes, respectively, where the different uncertain- 
ties obtained for the two decay modes are due to the 
different mass-window requirements. The uncer- 
tainty related with the w mass-window requirement 
is determined with the control sample J/q — wn 
with w — nnan’ and n — «*« «9. The dif- 
ference in w selection efficiency is estimated to be 
1.5% for the n. + wd mode. 


7. Background 


10. 


In the analysis of J/w — yod, the uncertain- 
ty associated with the peaking background from 
J/y > me, Ne > óK*K-, and 2(KtK7—) as 
well as the other unknown background is estimat- 
ed by varying up or down the numbers of back- 
ground events by one standard deviation accord- 
ing to the uncertainties of branching fractions in 
PDG B. The largest change in the ne > ġġ signal 
yield is determined to be 0.9%, and is taken as the 
systematic uncertainty. 


In the study of J/1» — ^w, the uncertainty asso- 
ciated with the unknown background is estimated 
by replacing the second-order Chebychev function 
with the first-order one. The change of the up- 
per limit of signal events is negligible. The uncer- 
tainty associated with the dominant background, 
Jl — nó — ^w, is estimated by varying the 
branching fraction by one standard deviation when 
normalizing the background in the fit. The differ- 
ence in the resulting upper limit is determined to 
be 5.696, and is taken as the systematic uncertainty. 


. Fit range 


In the nominal fit, the fit range is set to be Mgg 
and M. > 2.70 GeV/c?. Its uncertainty is esti- 
mated by setting the range of Mgg and Mug > 
2.60, 2.65, 2.75 or 2.80 GeV/c?. The branching 
fraction of ne — ġġ and the upper limit for ne > wọ 
are reestimated. The largest deviations to the nom- 
inal results, 0.7% for the decay ne — ġġ and 0.296 
for the decay ne — wọ, are taken as the systematic 
uncertainties. 


. Ne mass and width 


Uncertainties associated with the 7. mass and 
width are estimated by the alternative fits with the 
PDG values for the rj, parameters l. The resulting 
differences in the ne signal yield, 1.3% for n. > $$ 
and 5.6% for ne — wọ, are taken as systematic un- 
certainties. 


Amplitude analysis 

Systematic uncertainties associated with the am- 
plitude analysis arise including the uncertainties of 
the non-7- component and the mass resolution of 
Ne. 

In the nominal fit, the non-7. component is de- 
scribed by the nonresonant ¢¢-system assigned 
with quantum number J? = 07, 0+ and 2+. The 
statistical significance for the component with dif- 
ferent J? is determined according to the difference 
of log-likelihood value between the cases with and 
without this component included in the fit, taking 
into account the change in the number of degrees 
of freedom. The significances for the non-r com- 
ponent with J? = 07, 2+, 0* are 2.80, 3.00 and 
0.1c, respectively. If the 07 component is removed, 
the uncertainty is estimated to be +6.7%. If the 2* 
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TABLE I: Summary of all systematic uncertainties from the 
different resources (%). The combined uncertainty excludes 
the uncertainty associated with Dr(J/qp — yne), which is 
given separately. 


Sources Tl > 66 Ne — wó 
Nop 0.6 0.6 
Photon 1.0 3.0 
Tracking 4.0 4.0 
PID a 4.0 
Br(ó => K* K^) 2.0 1.0 
Br(w > nrn 79) — 0.8 
Kinematic fit 6.7 2.4 
Mia g= mass 0.7 L4 
M, +, zo mass m 1.5 
Background 0.9 5.6 
Fit range 0.7 0.2 
Ne mass and width 1.3 5.6 
Amplitude analysis Pii — 
Combined 1 10.7 
Br(Jfh — yne) 23.5 23.5 


component is removed, the uncertainty is estimated 
to be —26.096 mainly due to the strong interference 
between the n: and the 07 components. 


The uncertainty related with the 7. mass resolution 
is estimated by the alternative amplitude analysis 
with the detected width of the ne set to 34.2 MeV, 
estimated from the MC simulation with the nom- 
inal input ne width 32.0 MeV from Ref. [28]. The 
resulting difference of the ne signal yield with re- 
spect to the nominal value is 2.2%. 

'The total uncertainty from the amplitude analysis 
is estimated to be us 


Table |l] summarizes all sources of systematic uncer- 
tainties. The combined uncertainty is the quadratic 
sum of all uncertainties except for that associated with 
Br(J/b > Nc). 


VI. BRANCHING FRACTIONS 
A. ne> oo 


The product branching fraction of J/Y — yne — yoo 
is calculated by 


Na 
Br(J/v — yne)Br(ne — $6) = NP — K*K-) 


= (4.3 + 0.5(stat) "S (syst)) x 107°, 


where Br(ó — K*K-) is the branching fraction of 
the ç$ > K*K- decay taken from the PDG [2], Meg 
is the 7. signal yield, and e = 24% is the detection 
efficiency, determined with the MC sample generated 
with the amplitude model with parameters fixed accord- 
ing to the fit results. The number of J/wv events is 
Nj, = 223.7 x 109 [16]. 


9 


Using Br(J/ > yne) = (1.72:0.4)96 [2], Br(ne — $9) 
is calculated to be 


Br(n. = bo) - 
(2.5 + 0.3(stat) "O3 (syst) + 0.6(Br)) x 10 2, 


where the third uncertainty, which is dominant, is from 
the uncertainty of Br(J/w — Yne), and the second un- 
certainty is the quadratic sum of all other systematic un- 
certainties. 


B. m. — Gó 


No significant signal is observed for n. + wo, and we 
determine the upper limit at the 9096 C.L. for its branch- 
ing fraction, 


where N.p = 18 is the upper limit on the number of r) 
events at the 90% C.L., e = 5.9% is the detection efficien- 
Cy, Gsys = 25.8% is the total systematic error, and Br is 
the product branching fractions for the decay J/q) — Yne, 
$— K*K- and w — ata n? 


VII. SUMMARY AND DISCUSSION 


Using 223.7 million J/7 events accumulated with the 
BESIII detector, we perform an improved measurement 
on the decay of ne — ġo. The measured branching 
fraction is listed in Table and compared with the 
previous measurements. Within one standard devia- 
tion, our result is consistent with the previous measure- 
ments, but the precision is improved. No significant sig- 
nal for ne — wọ is observed. The upper limit at the 
9096 C.L. on the branching fraction is determined to be 
Br(ne — wo) < 2.5x 1075, which is 1 order in magnitude 
more stringent than the previous upper limit B. 

The measured branching fractions of ne + ¢¢ are three 
times larger than that calculated by next-to-leading per- 
turbative QCD (pQCD) together with higher twist con- 
hen ry This discrepancy between data and the 
HSR expectation If implies that nonperturbative mech- 
anisms play an important role in charmonium decay. To 
understand the HSR violation mechanism, a comparison 
between the experimental measurements and the theo- 
retical predictions based on the light quark mass correc- 
tion [13] the ? Po quark pair creation mechanism and 
the intermediate meson loop effects is presented in 
Table[IT] We note that the measured Br(n. — $49) is close 
to the predictions of the °P) quark model and the 
meson loop effects [14]. In addition, the measured upper 
limit for Br(n« + wọ) is comparable with the predicted 
value 3.25 x 10~4 in Ref. ug. The consistency between 
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TABLE II: Comparison of BESIII measured Br(n. — $$) with the previous results and theoretical predictions, where the 
branching fractions of n. — $$ from BESII and DM2 are recalculated with Br(J/4» — yne) = (1.7 + 0.4)96 | 


Experiment Br(J/v > ync)Br(ne > $$)(x10 ^) Br(n. — od) (x10 7) 
2.5 J d 


BESIII 4.3 £ 0.5703 E0.3793 £06 
BESII [5] 3.3 + 0.8 1.9 + 0.6 
DM2 [30] 3.9 + 1.1 2.3 + 0.8 
Theoretical Prediction Br(ne > $$) (x10 ^) 
PQCD[I] (0.7 ~ 0.8) 
? Py quark model [13] (1.9 ~ 2.0) 
Charm meson loop [14] 2.0 


data and the theoretical calculation indicates the impor- 
tance of QCD higher twist contributions or the presence 
of a non-pQCD mechanism. 
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Appendix 
Amplitude analysis of the decays 
Jf — yoo 


A. AMPLITUDES 


For the decay J/v(Ao) + Y(Ay)Ne > v60u)602), 
where the A;(i = y, 0, 1, 2) indicates helicity values for the 
corresponding particles, the helicity-coupling amplitude 
is given by 


An. (A0; A, M; À2) = FY (ri) DX: X, (00, G0) BW; (meg) 


i À F (Ey) 
x ae (r2) 95, X, (61, Qı) FE) , 
Yi 
(8) 


where rı(r2) is the momentum difference between y and 
Ne (two $ mesons) in the rest frame of J/w (ne), and 
Oo (ġo) and 04 ($1) are the polar (azimuthal) angles of 
the momentum vectors P4, and P, in the helicity sys- 
tem of J/w and ne, respectively. The z-axis defined for 
Ne > $()6(A2) is taken along the outgoing direction of 
$(A1) in the ne rest frame, and the z-axis is in the P}, 
and Pga) plane, which together with the new y-axis 
forms a right-hand system. BW;(m) denotes the Breit- 
Wigner parametrization for the 7. peak. The damp- 
ing factor F(E,) is taken as F(E,) = exp(— a) with 
8—0.065 GeV [26]; E? is the photon energy correspond- 


ing to Mog = mo. The helicity-coupling amplitudes F i 


and FY° ,, are related to the covariant amplitudes in the 
LS-coupling scheme by 


Bi(r1) 
py = op? = Sl, 1 
i Aa ` Bi) 
gu. Bilra) 
Fire = =p _ = T , 
i "U 2 Bi(r3) 
Foo = 0, (9) 


where B,(r) is the Blatt-Weisskopf factor [25], r? and 
r9 indicate the momentum differences for the two decays 
with mgs = Mny., and gis and gj, are the coupling con- 
stants for the two decays. 
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For the direct decay J/w — yoo, the mass spectrum 
of ġġ appears as a smooth distribution within the 7, sig- 
nal region; hence the Breit- Wigner function is excluded. 
The amplitudes for the direct decay are decomposed into 
partial waves associated with the ¢¢-system with quan- 
tum numbers J? = 07, 0* and 2*, and the high spin 
waves are neglected. These amplitudes are taken as 


Alyr(Ao, Ay, A1, A2) = FR, o Dyo,- (80; d0) FX, a; 
x D, as (61,01) for 07, 
ARA as d2) = FË o DS a, (80; do) FR a; 
ài Di X. (63, 01) for OF, 


F * 
AN R(A0, A; A ,À2) = 5 VFR. y DX xr, (90; 0) 
Àj 


nd * 
x FR DX ca (01, $1) for Qt. 


Here, helicity-coupling amplitudes F' FA are related to 
covariant amplitudes. For J? = 07, helicity amplitudes 
take the same form as that in Eq. (9). 

For the 0* case, helicity amplitudes are taken as 


FY = FY = A 
DUCTU EU A 
£ 2 B / 
[20 2 pot = 92272 2(r2) + Joo (10) 


+ 2 Do(r2) _ 90 
0 2 2V2 00 
Foo = "ES Ba(r9) n 3 


Fb = F°, 25 + 


Fü = pe e ON 0 
v35 Bar?) v5 Balr?) 
Bo(ri) , 9m 
2 D2\'1 0 
E Vas?" mo VIO 
(12) 


3 Ba(ri) gar? Bo(r1) 
Fo = F“, = V ag 94371 oy + - 0 


+ P j B»(ri) + 2923r? Bo(ri) | gor 
z 22 € 
27" | Bo(r) 35 Ba(r]) 30 
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3 Ba(r) gbor2 Ba(r2) 
F? = pe = UM ud pt 4 22072 
11 1-1 3524? Ba(r’) V3 Ba(r9) 
_ 93273 B>(r2) 4 902 
V21 Də(r2) v30 
(14) 
2 B (r) 1 Bo(r2) 
pot _ p? = 1 4 223 ea 2 
10 —10 vo Bir) 2921 2 Ba (r°) 
9273 Balr2) | 902 
2/7 Balr3) viO 
(15) 
^ Ba(r) 2 Bo(r2) 
p E p — 8427 4 Pa =; Pe 2iT2 
1-1 11 70 Bir!) 7922 2 Bx(r9) 
1 
Je 202. (16) 


For these nonresonant decays, the differences of momenta 
r) are calculated at the value m¿¿ = 2.55 GeV. 


The total amplitude is expressed by: 
A(Xo, Ay, A1;A2) = An. (A0, Ay, À1, A2) 


+ O A ayaa), (17) 
JP 


where the sum runs over J^ = 07,0* and 2*, and the 
symmetry of identical particle for two ¢ mesons is im- 
plied by exchanging their helicities and momentum. The 
differential cross section is given by 


3 
w= (5) 3. AQo, Ay, A1, Aa) 


Ña A EA AS 
x A*(Ao, Ay, À1, A2)d®, (18) 


where Ag, Ay = £1, and 1, à2 = +1,0, and d@ is the 
element of standard three-body phase space. 


B. FIT METHOD 


The relative magnitudes and phases for coupling con- 
stants are determined by an unbinned maximum likeli- 
hood fit. The joint probability density for observing N 
events in the data sample is 


N 
L= lI pa), (19) 


where P(z;) is a probability to produce event i with a set 
of four-vector momentum z; = (py, po, po)i. The normal- 
ized P(z;) is calculated from the differential cross section 
(dU /d9); 
TMC 


P(z;) = (20) 
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where the normalization factor omc is calculated from a 
MC sample with Nmc accepted events, which are gener- 
ated with a phase space model and then subject to the 
detector simulation, and are passed through the same 
event selection criteria as applied to the data analysis. 
With a MC sample of sufficiently large size, the cuc is 


evaluated with 
Nuc 
1 dr 
— — 21 
Nuc > s), QU 


OMC = 
For technical reasons, rather than maximizing £, S = 
— In Z is minimized using the package MINUIT. To sub- 
tract the background events, the In £ function is replaced 
with 


In £ = In Zaata — In Log. (22) 


After the parameters are determined in the fit, the 
signal yields of a given resonance can be estimated by 
scaling its cross section ratio R; to the number of net 
events, t.e., 


N; = R; (Nobs — Nyg), with R; — B $ (23) 
Ito 
where T; is the cross section for the ¿th resonance, Pot is 
the total cross section, and Nops and Npg are the numbers 
of observed events and background events, respectively. 
The statistical error, ó N;, associated with signal yields 
N, is estimated based on the covariance matrix, V, ob- 
tained from the fit according to 


Npars Npars 


ON; ON; 
P ct i i 
S X (mea), uec 2) 


m=1 n=1 


where X is a vector containing parameters, and u con- 
tains the fitted values for all parameters. The sum runs 
over all Npars parameters. 


C. RESULTS OF PARAMETERS 


The nominal fit includes the decays, J/ — Yne > 
yoo and J/v > «(6);» — yoo with JP = 07,0*,2*. 
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